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ORIGINAL ARTICLE 

Transcriptional regulation of the MET receptor tyrosine kinase 
gene by MeCP2 and sex-specific expression in autism and 
Rett syndrome 

JT Plummer 1 , OV Evgrafov 1 ' 2 , MY Bergman 3 , M Friez 4 , CA Haiman 5 , P Levitt 1 ' 6 and KA Aldinger 1 ' 7 

Single nucleotide variants (SNV) in the gene encoding the MET receptor tyrosine kinase have been associated with an increased 
risk for autism spectrum disorders (ASD). The MET promoter SNV rs1 858830 C low activity' allele is enriched in ASD, associated 
with reduced protein expression, and impacts functional and structural circuit connectivity in humans. To gain insight into 
the transcriptional regulation of MET on ASD-risk etiology, we examined an interaction between the methyl CpG-binding 
protein 2 (MeCP2) and the MET 5' promoter region. Mutations in MeCP2 cause Rett syndrome (RTT), a predominantly female 
neurodevelopmental disorder sharing some ASD clinical symptoms. MeCP2 binds to a region of the MET promoter containing the 
ASD-risk SNV, and displays rs1 858830 genotype-specific binding in human neural progenitor cells derived from the olfactory 
neuroepithelium. MeCP2 binding enhances MET expression in the presence of the rs1 858830 C allele, but MET transcription is 
attenuated by RTT-specific mutations in MeCP2. In the postmortem temporal cortex, a region normally enriched in MET, gene 
expression is reduced dramatically in females with RTT, although not due to enrichment of the rs1 858830 C low activity' allele. 
We newly identified a sex-based reduction in MET expression, with male ASD cases, but not female ASD cases compared with 
sex-matched controls. The experimental data reveal a prominent allele-specific regulation of MET transcription by MeCP2. 
The mechanisms underlying the pronounced reduction of MET in ASD and RTT temporal cortex are distinct and likely related 
to factors unique to each disorder, including a noted sex bias. 
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INTRODUCTION 

Autism spectrum disorders (ASD) are common, heterogeneous 
neurodevelopmental disorders that are clinically defined by 
impairments in social behavior, communication, and restricted 
interests and repetitive behaviors. 1-3 Common and rare single 
nucleotide variants (SNV) in many genes have been implicated as 
potential risk factors for ASD. 4-9 Demonstrating ASD causality due 
to coding mutations can be difficult, but variants in regulatory 
regions of risk genes pose additional challenges to investigating 
the neurobiological consequences of and contribution to ASD 
pathophysiology. The MET receptor tyrosine kinase is one gene for 
which an ASD-associated SNV has been shown to have functional 
consequences for gene function. A common SNV in the MET 5' 
promoter region (rs1 858830) was originally associated with ASD 
among families with multiple children diagnosed with ASD. 10 This 
finding has since been replicated in additional case-control and 
family association studies. 11-13 Additionally, MET expression is 
reduced in the temporal lobe of subjects with ASD. 14,15 More 
recently, neuroimaging studies demonstrated that the rs1 858830 
C allele impacts functional activity and structural connectivity in 
regions involved in social cognition in typically developing 
subjects, with a more pronounced effect in individuals with 



ASD. 16 The rsl 858830 C allele has also been associated with 
reduced gray matter growth in typically developing children and 
adolescents, 17 consistent with findings in mouse models that 
reducing MET expression disrupts neuronal architecture 18 and 
functional connectivity. 19 These data converge on the original 
finding that the rsl 858830 C allele reduces both nuclear protein 
binding to the MET promoter, and transcriptional activation of 
MET. 10 Furthermore, brain and peripheral MET protein levels are 
significantly lower in the presence of the rsl 858830 ASD-risk C 
allele compared with the non-risk G allele. 14,20 

Like other receptor tyrosine kinases, transcriptional regulation of 
MET expression is important for both normal and disease 
processes. 21,22 A number of transcription factors and DNA 
methylation patterns have been attributed to regulation of MET in 
cancer 22-25 but little is known regarding MET regulation in neural- 
relevant contexts. MET is a part of the biological network that 
includes several ASD-associated transcriptional regulators, including 
FOXP2 and MeCP2 26 FOXP2 mutations increase risk for language 
disorders 27 and direct FOXP2 binding to the 5' regulatory region of 
MET represses MET transcription. 28 Alterations in MECP2 cause severe 
neurodevelopmental disorders including Rett syndrome (RTT) and 
rare cases of ASD. 29-31 Given these findings, we hypothesized that 
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MeCP2 was a strong candidate as an additional transcriptional 
regulator of MET. The present report provides multiple lines of 
evidence that support this hypothesis, and further uncovered a 
previously unrecognized ASD sex-based and RTT-associated 
disruption of MET expression in human neocortex. 



MATERIALS AND METHODS 

Cultured olfactory neuroepithelial cells 

Protocols were approved by the Institutional Review Board at the University 
of Southern California and written informed consent was obtained from 
each subject. Genomic DNA samples from 27 control individuals that also 
had nasal biopsy tissue samples collected and olfactory neuroepithelial 
cultures established 32 were genotyped for rs1 858830 as described below. A 
total of nine cultured olfactory neuroepithelial cells (CNON) cells from male 
participants with mixed ancestry and representative rs1 858830 genotypes 
were cultured as previously described. 32 

Chromatin immunoprecipitation 

Human embryonic kidney (HEK) cells and CNON cells {n = 9) were grown to 
a confluence of 1 x 10 7 cells on 10 cm dishes as previously described 32,33 
Chromatin immunoprecipitation (ChIP) assays were conducted as 
previously described for HEK cells. 34 Approximately 25% of cells were 
pelleted and frozen for later RNA extraction. ChIP assays using CNON cells 
differed only in sonication time. Lysates from ON cells were sheared by 
sonication for a total of 8 min with pulsed intervals of 1 5 s ON and 45 s OFF 
on ice. For qPCR followed by ChIP, the LightCycler FastStart DNA Master PLUS 
Kit (Promega, Madison, Wl, USA) was used. Primers spanning 1.43 kb of 5' 
MET promoter were used for the ChIP qPCR (Supplementary Table S1). 
Assays were analyzed using a CFX96 real-time PCR detection system (Bio- 
Rad, Hercules, CA, USA). Fold enrichment was calculated relative to the 
immunoglobulin G ChIP and percent recovery was calculated relative to 
sample input. Total RNA was isolated from frozen CNON cell pellets (n = 6) 
and semi-quantitative real-time PCR was performed as previously 34 
described using rs1 858830 genotyping primers. 10 

Plasmid constructs 

Luciferase reporter plasmids pGL4.10 (empty) and pGL4.10[luc2] contain- 
ing 0.66 kb of MET promoter were previously described. 10 Coexpression of 
MeCP2 was accomplished using a MeCP2 cDNA clone (HsCD00081627) 
purchased from the DNASU repository (Biodesign Institute, Arizona State 
University, Tempe, AZ, USA). PCR was used to generate site-specific 
mutations in MeCP2 cDNA. PCR mutagenesis was performed according to 
the method described in the QuikChange Site-Directed Mutagenesis Kit 
(Stratagene, La Jolla, CA, USA). We utilized the MeCP2 cDNA plasmid as 
template and followed the manufacturer's primer design software 
(Stratagene). Primers are provided in Supplementary Table SI. All PCR 
was performed using Pfu Turbo (Stratagene) by initially denaturing the 
template at 95 °C for 30 s, followed by denaturing at 95 °C for 30 s, 
annealing at 60 °C for 1 min, extension at 68 °C for 7 min, with this cycle 
repeated 18 times. Original template DNA was digested by Dpn I 
treatment at 37 °C for 2h. Digested DNA was transformed into XL1-Blue 
cells for blue-white screening (Stratagene). Positive clones were purified 
using Promega Wizard Purification Kit (Promega, Madison, Wl, USA). The 
expected MeCP2 mutations were verified by DNA sequencing. 

MET 5' promoter luciferase assays 

HEK cells were plated onto 12-well plates. Twenty-four hours post plating, 
4(il of Lipofectamine2000 (Invitrogen, Carlsbad, CA, USA) was added to 
600 |il DMEM media containing a total of 0.3 [ig per well of the desired 
pGL4.10 reporter construct and/or MECP2 expression construct and 0.5 |ig 
of the reference Renilla luciferase reporter was included. The Lipofecta- 
mine and DNA solutions were then combined following manufacturer 
recommendations. Following 24 h of culture, cell lysates were prepared 
according to the manufacturer's recommendations of the Dual Luciferase 
Reporter Assay System (Promega). Both firefly and Renilla luciferase 
products were measured in the Tecan Infinite 200 Plate Reader (San Jose, 
CA, USA). Firefly luciferase activity was compared with Renilla luciferase 
activity as a relative ratio. This ratio represents the transcriptional activity of 
a particular luciferase reporter construct. Fold of activation in the luciferase 
assays was calculated after normalization against the empty firefly control 
vector. Each experiment was performed minimally in triplicate. 



RNA and DNA isolation from human postmortem brain samples 
Fresh-frozen postmortem brain samples were obtained through the 
Autism Speaks-supported Autism Tissue Program at the Harvard Brain 
Tissue Resource Center (http://www.brainbank.mclean.org/) or the NICHD 
Brain and Tissue Bank at the University of Maryland School of Medicine 
(http://www.medschool.umaryland.edu/BTBank/). Superior temporal gyrus 
samples were obtained from 15 ASD and 5 RTT brain samples that were 
sex, age and postmortem interval (when possible) matched to 18 CTL brain 
samples. 34,35 The majority of samples (37/38) were from individuals of 
European descent; one sample was from an individual of African-American 
descent. Total RNA was isolated using the mirVana miRNA Isolation Kit 
(Invitrogen) according to the manufacturer's protocol. Total RNA 
concentration was determined using an Agilent Bioanalyzer 2100 system 
(Santa Clara, CA, USA); average RNA integrity number (RIN) = 7.3. Genomic 
DNA isolation was performed previously. 34,35 

MET rs1 858830 genotyping in RTT and control samples 
Genomic DNA was obtained for 193 unrelated females positive for MECP2 
mutation or exonic deletion from the Greenwood Genetic Center. An 
additional 15 genomic DNA samples from RTT females were purchased 
from the Coriell Cell Repository (http://www.ccr.coriell.org/). Genomic DNA 
samples from 514 unrelated females from the Multiethnic Cohort were 
used as controls 36 DNA was amplified with the KOD Xtreme Hot Start PCR 
kit (EMD Biosciences, San Diego, CA, USA) according to the manufacturer's 
protocol using published primers. 10 Amplicons were submitted for direct 
resequencing (BeckmanCoulter Genomics, Beverly, MA, USA or Eton 
Biosciences). Sequence electropherograms were analyzed using 
Sequencher v5.0 (Gene Codes Corporation, Ann Arbor, Ml, USA). Autism 
Genetic Resource Exchange (AGRE) rsl 858830 genotypes were previously 
determined 10,11 and downloaded with permission from the AGRE website 
(https://research.agre.org). Population data for rsl 858830 were downloaded 
from the 1000 Genomes Browser (http://browser.1000genomes.org 
(August 2012)). Allele and genotype data were analyzed using % 2 and a 
two-sided P<0.05 was considered significant. Allelic % 2 results are 
presented. The Genetic Power Calculator (http://pngu.mgh.harvard.edu/ 
~purcell/gpc/) 37 was used to estimate the power of the case-control 
association analyses. 

Quantitative real-time PCR postmortem brain 
cDNA was synthesized from total RNA using the Superscript III First Strand 
cDNA Synthesis Kit (Invitrogen) according to the manufacturer's protocol 
using oligo dT primers. TaqMan Gene Expression assays for MET (assay ID 
Hs01565584_m1), MECP2 (assay IDs HsOOl 72845_m1 (M£CP2_e1) and 
Hs01 598237 {MECP2_e2)), and glyceraldehyde-3-phosphate dehydrogen- 
ase(assay ID Hs9999905_m1) were purchased from Life Technologies 
(Carlsbad, CA, USA). qPCR was performed in triplicate for each sample 
using lOOng cDNA and TaqMan Gene Expression Master Mix (Applied 
Biosystems, Foster City, CA, USA) according to the manufacturer's protocol. 
Assays were analyzed using a CFX Real-Time PCR detection system (Bio- 
Rad). The delta cycle threshold (AC t ) of target relative to glyceraldehyde-3- 
phosphate dehydrogenase was averaged for each sample. Analysis of 
variance was performed using AC t values. Results are presented as relative 
expression of target compared with glyceraldehyde-3-phosphate 
dehydrogenase. 

RESULTS 

MeCP2 binds to the MET 5' promoter ASD-associated region 
We first identified MeCP2 as a putative regulator of MET 
transcription during an in vitro TF screen using arrays containing 
140 TF proteins (Supplementary Figure S1). An oligonucleotide 
probe generated from the 5' region of MET bound to several TF 
proteins, including Hand2, Lhx2 and MeCP2. We further investi- 
gated the putative binding of MeCP2 to the MET promoter using 
several additional assays. MeCP2 binds to CpG dinucleotides, 
having a complex role in transcriptional regulation during brain 
development. 38,39 The MET 5' promoter harbors a putative 700 bp 
CpG island containing >70 CpG sites (Supplementary Methods 
and Figure 1a). We examined MeCP2 binding within the 5' MET 
promoter in HEK cells by ChIP. Anti-MeCP2 specifically precipitates 
MET DNA, whereas control immunoglobulin G yields no MET 
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Figure 1. MeCP2 directly binds to the MET promoter, (a) Schematic of the 5' promoter region of MET drawn to scale; Human Genome Browser 
(hg19), chr7: 116098419-116099867. Horizontal lines indicate the relative locations for the CpG island (green), primers used in ChIP assays 
(gray), functional promoter variant (rs1 858830), MET transcriptional start site (TSS). (b) Anti-MeCP2 antibody directly pulls down the MET 
promoter sequence using primers 2 and 7 (red bars). Anti-acetyl histone 3 (H3) and anti-immunoglobulin G (IgG) antibodies were used as 
positive and negative controls, respectively, (c) qPCR analysis of ChIP by MeCP2 of the MET promoter sequence using primer 7. *P< 0.001. 



enrichment. PCR using primers tiled across the MET 5' promoter 
region revealed that MeCP2 binds to multiple regions (Figures 1a 
and b). By ChIP, we show robust MeCP2 binding within the CpG 
island (primers 2 and 7). Primer 7 includes the region that contains 
the rsl 858830 ASD-risk SNV (Figures 1a and c and Supplementary 
Figure S2). 

The rsl 858830 C allele enhances MeCP2 transcriptional regulation 
of MET 

To determine whether the rsl 858830 SNV influences MeCP2 
transcriptional regulation of MET, we transfected two luciferase 
reporter constructs containing 663 bp of the MET promoter 
(Figure 1a), differing only at the rsl 858830 nucleotide, together 
with MECP2 cDNA into HEK cells. First, we replicated previous 
findings 10 that the reporter construct containing the G allele 
generates greater luciferase activity compared with the construct 
containing the C allele (P = 0.022; Figure 2a). Next, cells were 
cotransfected with MECP2 cDNA and each of the two MET 
promoter constructs. Surprisingly, greater luciferase activity was 
detected when the C allele construct was cotransfected with 
MECP2 compared with the C allele alone (P = 0.0002). No 
significant difference in luciferase activity was detected when 
the G allele was cotransfected with MECP2 compared with the G 
allele alone (P= 0.099). In the presence of MeCP2, the C allele also 
showed greater transcriptional activity than the G allele 
(P = 0.0001). These data indicate that the rsl 858830 C allele can 
directly modulate MeCP2 activation of MET transcription. 



mutations could disrupt MET promoter activity. Several constructs 
containing common RTT-causing MECP2 mutations located within 
the methyl-binding domain or the transcription repressor domain 
of MeCP2 (Figure 2b) were generated using site-directed 
mutagenesis. The mutant constructs were transfected into HEK 
cells and luciferase activity was monitored to assess transcriptional 
activity of the MET promoter (Figure 2c). Again, cotransfection of 
MECP2 cDNA and the rsl 858830 C allele construct activated MET 
transcription. A two-way ANOVA followed by Tukey HSD post hoc 
comparisons were used to evaluate the statistical significance 
of rsl 858830 allele by MECP2 mutation on MET transcription. 
Both rsl 858830 allele and MECP2 mutation had a significant 
effect on MET transcription {P = 0.038 and P = 0.002, respectively). 
The interaction between rsl 858830 allele and MECP2 mutation 
approached significance {P= 0.055). In comparison with wild-type 
MECP2, MECP2 with p.T158M or p.R168X mutation failed to 
enhance the activation of the C allele (P = 0.017 and P = 0.048, 
respectively). Notably, no mutations in the methyl-binding domain 
of MeCP2 significantly altered MET transcription in the presence of 
the G allele compared with wild-type MECP2 (P^ 0.879). Mutations 
in the transcription repressor domain of MeCP2 (p.R270X and 
p.R306C) did not significantly alter MET transcriptional activation 
in the presence of either the C or G allele (P^ 0.627). Thus, in the 
presence of the rsl 858830 C allele, MET transcription is attenuated 
by RTT-specific mutations in MeCP2 (p.T158M and p.R168X) that 
impact the methyl-binding domain, but not those in the 
transcription repressor domain. 



Overexpression of mutant MECP2 impacts MET transcriptional 
regulation 

To further establish a role for MeCP2 in the transcriptional 
activation of MET, we next tested whether RTT-causing MECP2 



MeCP2-MET binding is ASD-risk rsl 858830 genotype-dependent 
in neural progenitors 

To examine allele-specific binding at rsl 858830 by MeCP2, we 
performed ChIP on DNA isolated from cultured primary neural 
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Figure 2. Functional characterization of MeCP2 and MECP2 
mutations on MET transcriptional activation, (a) Luciferase reporter 
assays demonstrate differential activation of the MET promoter 
by MeCP2. MET luciferase reporter constructs containing rs1 858830 
G or C were transiently transfected into HEK cells with or 
without addition of MECP2 cDNA. *P< 0.001. (b) Schematic of the 
MeCP2 protein structure with common RTT-causing mutations 
(MBD; TRD, transcription repressor domain; NLS, nuclear localization 
signal), (c) Luciferase assays of MECP2 mutations cotransfected with 
MET promoter luciferase showed altered transcription compared 
with wild-type MECP2. (red line-rsl 858830 C allele + WT MECP2; 
blue line-rsl 858830 G allele + WT MECP2). *P<0.05 compared 
with rs1 858830 allele (G or C respectively) + WT MECP2. Fold 
of luciferase activation was calculated after normalization 
against an empty luciferase control vector. All transfections 
were performed in triplicate. 



progenitor cells derived from human olfactory neuroepithelium of 
normal subjects followed by qPCR using primers within the region 
containing rs1 858830, as in the HEK ChIP experiments (Figure 3). 
Genomic DNA was previously isolated from peripheral blood 
samples from 27 control subjects who underwent a nasal biopsy. 32 
We determined the genotype at rs1 858830 by direct sequencing 
and selected three different CNON for each of the three rsl 858830 
genotypes. Anti-MeCP2 specifically precipitates MET DNA in the 
presence of the ASD-risk genotype (CC), but not in the presence of 
the ASD-non-risk genotypes (GG and GC) (P = 0.037). The negative 
control immunoglobulin G does not yield enrichment of the MET 
sequence (P = 0.642). MET expression was not significantly 
different among the ASD-risk (CC), the heterozygous (GC) 
genotype and the non-risk (GG) genotypes (P = 0.059; 
Supplementary Figure S3). However, we noted a trend for 
increased MET expression in the CNON cells with the rsl 858830 
ASD-risk (CC) genotype. 
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Figure 3. Differential binding of MeCP2 in the presence of rsl 858830 
ASD risk (C) and non-risk (G) alleles. Human CNON cells were 
assayed by ChIP using an anti-MeCP2 antibody. qPCR analysis using 
primer 7 showed increased MeCP2 binding to the MET promoter in 
CNON cells homozygous for the risk (CC) genotype compared with 
CNON cells heterozygous (GC) or homozygous for the non-risk (GG) 
allele. Results reflect the mean ± s.e.m. across three CNON cells for 
each genotype. P< 0.037. 



MET expression is reduced in postmortem brain tissue of females 
with RTT and males with ASD 

Few risk genes implicated in ASD have been directly examined in 
RTT cases. Thus, in order to translate the in vitro findings of a direct 
MeCP2-ME7 interaction, we assayed MET expression by qPCR in 
temporal cortex from age-matched RTT and control females 
(Figure 4a). Temporal cortex was used due to the enrichment of 
MET expression in this region of primate neocortex. 1 5 ' 28 ' 40,41 MET 
expression was reduced dramatically in the temporal cortex of 
females with RTT compared with sex-matched controls (P = 0.007). 

A significant reduction in MET protein was previously described 
in the temporal cortex of individuals with ASD, 14 with similar 
findings when transcript was examined; 15 these studies comprise 
predominantly male subjects. The RTT findings raised the 
possibility of differential sex-based changes of MET in ASD. Thus, 
we assayed MET expression by qPCR in temporal cortex from age- 
and sex-matched ASD and control individuals (Figure 4b). There 
was a significant reduction in MET expression in the temporal 
cortex of all individuals with ASD compared with controls 
(P = 0.032), consistent with previous results. Stratification by sex 
further showed a significant reduction in MET expression in males 
(P = 0.016), but not females (P = 0.720), with ASD compared with 
controls (Figure 4c). MET expression is not different between male 
and female controls (P = 0.403). 



Translational Psychiatry (2013), 1 -8 



© 2013 Macmillan Publishers Limited 



Role of MeCP2 in autism and Rett syndrome 
JT Plummer et al 



CTL 



ASD 



jjj 0.005 
a. 

& 

fc 0 004 
2 

■o 0.003 



E 0.002 
8 

z 0.001 
0 



o 0.006 

8 

£ 0.005 
| 

,!! o.oo4 
u 

^ 0.003 

s 

^ 0.002 

£ 0.001 
0 



8 


> 

ooo 


§ — 0 — 


n 


OO 


O 


O 


I* 


8 


— O — 

, B , 


CTL 


ASD 



§ 0.006 

05 

£ 0.005 
& 

£ 0.004 

Hi 

^ 0.003 
I 

| 0.002 

z 0.001 
0 



8 




O 


8 


— o — 


a 


Q 


— O — 


O 


o 


O 


O 


O 




o 






O 


O 


OC 


0 

o 


M 


F 


M 


F 



Figure 4. Sex-specific MET expression in postmortem brain of individuals with ASD. Samples are represented by open circles and group means 
are represented by horizontal bars, (a) MET expression in temporal cortex of females with RTT and controls (CTL). **P< 0.001 (b) MET 
expression in temporal cortex of individuals with ASD and controls, (c) MET expression in the temporal cortex of individuals with ASD and 
controls as shown in panel b, separated by sex. No significant difference in MET expression was detected between males (M) and females (F) 
among controls. *P<0.05. 



MET protein levels in the temporal cortex were previously 
associated with rs1 858830 genotype in controls. 14 Stratification by 
rs1 858830 genotype showed reduced MET expression in ASD 
across all genotypes compared with controls (Supplementary 
Figure S4). Although there was no significant effect of rs1 858830 
genotype on MET expression for ASD (P = 0.220) or controls 
(P = 0.969), the C allele was overrepresented among the brain 
samples (59.7%). After correcting for rs1 858830 genotype, MET 
expression remained significantly reduced in males (P = 0.021), 
and not females (P= 0.269), with ASD compared with controls. 
The C allele was also overrepresented in RTT. However, the lack of 
GG genotype precluded statistical power to perform the same 
analysis. 

No difference in MECP2 expression was detected in temporal 
cortex among cases and controls (Supplementary Figure S5), nor 
between males and females (Supplementary Figure S6). 

The ASD-risk rsl 858830 C allele is not associated with RTT 
Results from the in vitro luciferase experiments suggest MeCP2 
binds to the MET promoter regardless of which allele at the 
rsl 858830 nucleotide is present. However, when the rsl 858830 C 
allele is present, MeCP2 binding enhances MET transcription, 
unless the bound MeCP2 protein contains a mutation in the 
methyl-binding domain. These data predict MET should be 
expressed at normal levels in RTT, yet MET expression in temporal 
cortex is significantly reduced in RTT compared with controls. 
Previous studies have established that the rsl 858830 C allele 
decreases MET transcription and leads to reduced MET protein 
levels in postmortem brain and peripheral cells. 14,20 An association 
between MET expression and rsl 858830 could not be evaluated in 
RTT brain due to insufficient sample size. To address whether the 
low MET expression in RTT could be explained by an enrichment 
of the rsl 858830 C allele, we conducted a case-control study of 
females with RTT and unrelated, sex-matched unaffected subjects. 
Genotype at rsl 858830 was determined by direct resequencing in 
208 unrelated RTT and 514 control females. No significant 
differences in genotypic and allelic frequencies were detected 
between RTT and control females (Supplementary Table S2). The 
rsl 858830 C allele frequencies also were not significantly different 
in females with RTT compared with the mothers of children with 
ASD (P = 0.695), the parents of children with ASD (P = 0.786), nor 
the allele frequencies in the 1000 Genomes Project 42 (P = 0.176). 
The case-control analyses had 77-88% power to detect associa- 
tion in an allelic test with a relative risk of 2.27 for the CC genotype 
and 1.67 for the CG genotype, as reported in the initial association 
study. 10 Similar results were found in a European-only analysis 
(Supplementary Table S3), suggesting that the lack of association 



was not due to population stratification. Together, the expression 
and genotype data suggest that the mechanisms through which 
reduction of MET occurs in ASD and RTT are distinct and likely 
related to factors specific to each disorder. 



DISCUSSION 

The present study demonstrates binding of MeCP2 to the MET 
promoter that is functionally relevant, with regulation of MET 
transcription by MeCP2 influenced by the rsl 858830 SNV in vitro. 
This regulation is likely to be cell context specific, as binding 
differences were discovered between HEK cells and primary 
human neural progenitors from the olfactory epithelium. More- 
over, the analyses of postmortem brain samples from cases and 
controls for both ASD and RTT demonstrate sex-based differences 
in reduced MET expression, which has functional implications in 
light of the role of MET in cortical development and circuit 
function, 16,18,19 and in individuals with ASD. 43 Although our data 
also reveal pronounced reduction of MET in both ASD and RTT 
temporal cortex, the sex-based differences appear to occur 
through distinct mechanisms, with males impacted in ASD and 
females in RTT. Future analyses are required to provide additional 
insight with regard to mechanisms, but irrespective of these 
differences, these findings indicate that like other molecules 
involved in synaptic and circuit development, disruption of MET 
expression in different neurodevelopmental disorders, either 
directly through genetic variants that impact transcription (ASD), 
or indirectly through as yet unknown factors (RTT), will have 
functional consequences that contribute to disorder symptoms. 

Implications of MeCP2 Regulation of MET 

MeCP2 acts as a global transcriptional regulator by recruiting 
chromatin-remodeling complexes or regulating higher-order 
chromatin structures 44-50 Thus, MeCP2 function may be 
determined by its interaction with numerous protein partners 
that produce functional outcomes, or MeCP2 may globally alter 
chromatin state to regulate transcription based on the status of 
DNA methylation and MeCP2 activation. 51 Each of these scenarios 
is consistent with the modest changes in gene expression 
detected in human and mouse tissues with altered levels of 
MECP2. 52 ~ 60 These complex molecular and biochemical inter- 
actions obfuscate the impact of MeCP2 transcriptional regulation 
of specific genes. Thus, few specific MeCP2-regulated genes are 
known. 61-63 However, large number of genes are dysregulated in 
discrete brain regions of Mecp2 mouse models, 52,59,60 suggesting 
pronounced, yet restricted cell-specific influences of MeCP2 may 
depend on both developmental and physiological states. 50,64-66 
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These interactions would be difficult to discern in analyses of 
human brain tissue. 

The indirect connection between MET and MeCP2, first high- 
lighted in a network model of ASD-risk genes, 26 has been 
demonstrated as a direct relationship here using methodologies 
that measure both protein-to-DNA binding and transcriptional 
activity. MeCP2 binding to the MET promoter within CpG domains 
raises the possibility that there may be activity-dependent 
changes in DNA methylation status in combination with MeCP2- 
mediated chromatin state to regulate MET transcriptional 
regulatory complexes, as is seen with MeCP2 regulation in 
mouse brain development and function. 50,67,68 This is consistent 
with the recent discovery of environmental factors that alter MET 
protein expression in vitro and in vivo, 69,70 reflecting the regulatory 
sensitivity of this gene. The data measuring transcriptional activity 
of the MET promoter support this concept, as both the C and G 
allele are permissive for MeCP2-mediated transcriptional activity, 
but the magnitude of MeCP2 binding to the MET promoter was 
less pronounced in HEK cells relative to CNON cells. These findings 
indicate additional MeCP2 cofactors specific to neural progenitor 
cells are responsible for the differential activity of rs1 858830 G 
versus C. Identifying the specific protein complexes that are 
influenced by MeCP2 binding to the MET promoter are part of 
ongoing investigations. Additionally, differences in DNA methyl- 
ation at the MET promoter between HEK cells and CNON cells may 
contribute to the differential effects of MeCP2 transcriptional 
regulation detected here; future studies are required to test this 
possibility. 



Sex-Based Differences in MET Expression in Neurodevelopmental 
Disorders 

Perhaps most surprising in our studies was the finding that there 
are sex-based differences in MET expression in the temporal 
neocortex. This arose from our analysis of RTT postmortem 
neocortical tissue, which had not been investigated previously for 
MET expression. We based the rationale on two factors: (1) the 
cooccurrence of gastrointestinal disturbances in girls with RTT; 71 
the MET risk allele is enriched in children with ASD and gastro- 
intestinal disturbances; 72 and (2) in mice, Mecp2 and Met fore- 
brain expression peaks during periods of dendritic growth and 
synaptic formation, and each mouse mutant displays synaptic 
hyperconnectivity 65 reminiscent of MET influence on functional 
activation and network connectivity in humans 16 . Though there 
were a limited number of RTT postmortem brains for analysis, the 
nearly undetectable levels of MET compared with matched female 
controls was highly statistically significant and has functional 
implications. Because reports of altered gene and protein 
expression in postmortem cases of ASD are generally dominated 
by male subjects (due to the 5:1 ratio of male:female diagnoses), 
minimal data were available on MET expression in females with 
ASD. In contrast to what we found in RTT, MET expression was 
reduced in male, but not female temporal neocortex. There are 
several noteworthy conclusions from these data. First, RTT and 
ASD only share some common clinical symptoms during the 
initial regressive period of RTT 73 , suggesting there are distinct 
mechanisms that underlie the primary etiologies. Notably, a major 
difference between ASD and RTT etiology is deletion or mutation 
of MECP2 in RTT patients. However, RTT remains a clinical diag- 
nosis that displays a range of severity in its clinical presentation, 
suggesting additional factors modulate clinical characteristics. 
Our results suggest mutant MeCP2 binding fails to enhance MET 
expression in some cases, rather than decreases MET expression 
when the rs1 858830 C allele is present, yet MET expression 
is reduced in RTT brain. In ASD, MET expression is reduced 
in the presence of the rs1 858830 C allele 14,20 . Thus, the factors 
that influence the sex-based reductions of MET must be disorder 
specific. Second, the ASD-associated rs1 858830 C allele is 
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associated with severe social and communication phenotypes, 
and is enriched even further in children with ASD and co- 
occurring gastrointestinal disturbances, 72 yet it is not associated 
with RTT, which shows both severe communication and 
gastrointestinal disturbances phenotypes. There are several 
possible explanations. Other MET variants associated with 
ASD 13,74 or variants in genes that regulate MET transcription 
could be enriched in RTT cases. Alternatively, non-heritable factors 
may contribute to the pronounced reduction in MET. These 
analyses raise the issue of whether there are fundamental sex- 
based differences in normal transcriptional regulation of MET, as 
well as of other ASD and neurodevelopmental disorder risk genes. 
Ongoing studies examining larger human female and male 
cohorts, as well as experiments in model systems are required 
to determine sex-based differences of intrinsic regulation and 
response to environmental factors in order to elucidate novel 
mechanisms relevant to understanding disorder etiologies. 
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